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We consider a two-parameter one-dimensional Hamiltonian with uncorrelated diagonal disorder 
and non-random long-range inter-site interaction J mn = J/\m — . The model is critical at 
1 < fi < 3/2 and reveals the localization-delocalization transition with respect to the disorder 
magnitude. To detect the transition we analyze level and wave function statistics. It is demonstrated 
also that in the marginal case (fi = 3/2) all states are localized. 
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Localization-delocalization transition (LDT) in disor- 
dered systems, predicted by Anderson for three dimen- 
sions (3D) in 1958,» (see also Ref. 2) still remains a fas- 
cinating problem (see Refs.00H for an overview). Dur- 
ing the last two decades, a remarkable progress has been 
achieved in understanding the LDT, especially in discov- 
ering the nature of the wave function at transition. This 
progress became possible thanks to the fruitful idea of 
the multifractality of wave functions at criticalityi£i2'2iill 
This conjecture was then analytically proven for an en- 
semble of power-law random banded matrices (PRMB), 
which revealed the LDT with respect to the interaction 
exponen iLLLi (see Ref. for an overview). Within the 
framework of the latter, it was demonstrated, in par- 
ticular, that (i) the distribution function of the inverse 
participation ratio (IPR) is scale invariant at transition 
and (ii) the relative IPR fluctuation (the ratio standard- 
deviation/mean) is of the order of unity at the critical 
point This finding confirmed the conjecture, that 
was put forward for the first time in Refs. that 
distributions of relevant physical magnitudes are univer- 
sal at criticality (see also Refs. lr^lT^llai . This invariance 
is then a powerful tool to monitor the critical point. 

In the present paper, we consider a two-parameter 
tight-binding Hamiltonian on a regular ID lattice of size 
iV with non-random long-range inter-site interaction: 
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where |n) is the ket vector of a state with on-site energy 
e n . These energies are stochastic variables, uncorrelated 
for different sites and distributed uniformly around zero 
within the interval of width A. The hopping integrals 
are J mn = J/\m — n| M , J nn = with 1 < fi < 3/2. For 
definiteness we set J > 0, then the LDT with respect 
to disorder magnitude occurs at the upper band edge, 



provided 1 < fi < 3/2i2Si2i The transition is analogous 
to that within the standard 3D Anderson model, /i = 
3/2 represents the marginal case in which all states are 
expected to be weakly localized. 21 

To detect the transition we analyze level and wave 
function statistics. We perform a numerical analysis of 
size and disorder scaling of the relative fluctuation of 
both the nearest-level spacing (LS) and the participation 
number (PN). The latter is defined as: 
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where %j) vn denotes the n th component of the v th nor- 
malized eigenstate of the Hamiltonian 

The relative fluctuation of the nearest-level spacing is 
an invariant parameter at transition, as was conjectured 
in Ref. |u] for the 3D Anderson model and demonstrated 
later for a variety of other disordered models (see e.g. 
Refs. and references therein). The invariance can be 
used to detect the critical point. We demonstrate that 
within the present model, the ratio of the standard devi- 
ation of the PN (SDPN) to its mean value (MPN) is also 
an invariant parameter at the critical disorder magnitude 
A c . Therefore, the ratio SDPN/MPN can also be used 
to detect the transition. To the best of our knowledge, 
this quantity has never been used for this purpose. 

As the LDT occurs at the top of the band within the 
considered two-parameter model, we calculate disorder 
and size scaling for uppermost states. Open chains are 
used in all calculations. We take advantage of the Lanc- 
zos method to calculate the scaling for large system sizes 
(up to about 6 x 10 4 sites) and two particular values of 
the interaction exponent: \i = 4/3 (the LDT occurs) and 
ji = 3/2 (the marginal case; no transition is expectedSi). 
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FIG. 1: Disorder scaling of the relative fluctuation of the 
nearest-level spacing (the ratio SDLS/MLS) for /i = 4/3 in 
the vicinity of the joint intersection point (that is at A = 
10.7 -r 11.5 J). The curves are calculated for different system 
sizes TV and averaged over more than 5 • 10 3 x (65536/7V) 
disorder realizations. 
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FIG. 2: Disorder scaling of the relative fluctuation of the PN 
(the ratio SDPN/MPN) for n = 4/3 in the vicinity of the joint 
intersection point (that is at A = 10.0 -r 10.6 J). The curves 
are calculated for different system sizes N and averaged over 
more than 5 • 10 3 x (65536/7V) disorder realizations. 



First, we calculate the critical point by means of the 
level statistics analysis. In Fig. ^ we plotted the dis- 
order scaling of the ratio of the standard deviation of 
the nearest-level spacing (SDLS) distribution to its mean 
(MLS) at the top of the band for \i = 4/3. The figure 
demonstrates that all disorder-scaling curves plotted for 
different system sizes intersect within a narrow range of 
A, between 10.7 J and 11.5J. 

Calculations of the scaling of the relative PN fluctua- 
tion confirm the conjecture that the ratio SDPN/MPN is 
also a size invariant parameter at transition: Fig.[5]shows 
that all SDPN/MPN curves plotted versus disorder for 
different system sizes intersect in a narrow range of A, 
from 10.0 J to 10.6 J. One can deduce from Fig. [21 that 
both the MPN and the SDPN are of the same order of 
magnitude at the intersection fo r any system si ze, as was 
shown for other models in Refs. lllll3ll4ll8lT9l 

The regular size dependence of intersection points in 
Figs.n]and|2|is a finite size effect; accounting for the lat- 
ter by means of the finite size scaling analysis allows for 
obtaining the value of the critical disorder. Both figures 
demonstrate that finite size effects are unusually strong 
which, within the present model, results from the long- 
range nature of the inter-site interaction. Contrary to 
the standard Anderson model, the contribution of the 
long-range coupling terms to the spectrum of the Hamil- 
tonian converges very slowly as the system size in- 
creases. The latter results in a corresponding increase of 
the band width (mostly, the upper band edge, where the 
LDT takes place). For an open chain, the upper band 
edge E(N) size-scales as follows: 

E{N) = E 00 ( fJ ,)-^ I +0(N-n. (3) 



For fx = 4/3, £^(4/3) « 7.20 J and C(4/3) * 8.45 J. The 
increase of the band width with the system size leads to 
the fact that disorder of the same magnitude is effectively 
weaker for larger systems. The latter effect introduces 
regular size dependence of the critical disorder that is 
obtained by numerical analyses of finite systems. The 
contribution of other finite size effects^LS such as in- 
fluence of boundary regions, are expected to be weaker 
for large systems because of very slow convergence of the 
upper band edge (oc ./V 1_M ). Our calculations confirm 
this conjecture. 

The intersection point, A(Ni , N2), of two disorder scal- 
ing curves plotted for different system sizes N± and N 2 
depends on the sizes. To account for such dependencies 
we proceed as follows. First, set by definition: 

A C (N) = A(N-1,N+1), N r ^oo, (4) 

where A C (N) is the critical disorder that can be obtained 
by analyses of a finite system of size N (N r is the num- 
ber of disorder realizations over which the averaging is 
performed). Second, use the following anzats for the in- 
tersection point: 

A(JVi, N 2 ) = w (Ni.) A c (Nt) + w{N 2 ) A C {N 2 ) , (5) 

where the weight function w(N) is to be determined. 
Bearing in mind the slow convergence of the band edge 
(oc N 1 ^^), we use the following anzats for A C (N): 

A C (N) w A c (oo) + bN 1 ~^ + cN-~ t , 7V>1, (6) 

where the 6, c, and 7 > ji — 1 arc fitting parameters. 
Using the anzats © together with Eqs. (JSJ and (JSJ and 
expanding A(7V-1, N+l) in series about N (at N > 1), 
we find the weight function: w(N) = 1/2 + 0(1/A^ +1 ). 
Further, for any given pair Ni < N 2 (Ni, N 2 3> 1) there 
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exists the size N, such that A C (AT) = A(N l ,N 2 ). Mak- 
ing use of the latter equation together with Eq. iJBJ and 
keeping the leading (non-zero) power of system size in all 
expansions, we find the sought N: 
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Thus, the intersection point of disorder-scaling curves 
plotted for two different system sizes Ni < N 2 (Ni , N 2 
1) yields the critical disorder for an intermediate system 
size N as defined by Eq. JJJ). 

We further use Eq. (|7|) and intersection points of the 
curves in Fig. □ (LS data) and Fig. (PN data) to ob- 
tain A c (TV). Figure shows A C (N) together with best 
nonlinear fits of Eq. @ to the whole data sets (dashed 
lines) and the best linear fit of E(N) given by Eq. (J3J) to 
the three last PN-data points (solid line). The nonlin- 
ear fits give A c (oo) = (10.97 ± 0.09) J for the LS data, 
and A c (oo) = (11.19 ± 0.10)J for the whole PN data 
set, while the linear fit of E(N) to the tail PN points 
gives A c (oo) = (10.91 ± 0.17) J. The obtained values of 
A c (oo) agree well with each other. This confirms our 
conjecture that for large system sizes the band edge size 
dependence provides the dominant contribution to the 
finite size effects. Finally, the critical disorder is deter- 
mined as A c (oo) = (11.09 ± 0.21) J for fi = 4/3. 
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FIG. 3: Critical disorder sizes scaling obtained from the PN 
data (■) and level statistics data (k). Dashed lines are best 
fits of Eq. @ to the LS data (b = -26.44, c = 824.16, 7 = 
0.70) and the PN data (b = -21.61, c = 12.73, 7 = 2.82). 
The solid line is the best linear fit of Eq. © to the last three 
PN data points: A C (N) = (1.52 ± 0.02) x E(N). 

It should be noticed that, despite that both methods to 
detect the LDT are in good agreement, finite size effects 
are more pronounced in the case of analysis of the level 
statistics. For these reasons, the proposed method which 
is based on the studies of the wave function statistics 
appears to be advantageous, at least for the considered 
model. 
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FIG. 4: a) Disorder scaling of the relative PN fluctuation 
(the ratio SDPN/MPN) for fj, = 3/2 in the vicinity of the joint 
intersection point at A c = 0. The curves are calculated for 
two different system sizes (65536 and 8192) and averaged over 
more than 5 x 10 3 and 10 5 disorder realizations respectively, 
b) A blow up of the crossing at the origin. 



We applied the same technique to analyze the local- 
ization properties in the marginal case, fj, = 3/2, where 
the states are expected to be localized weaklyiSi Figure^] 
shows the SDPN/MPN scaling curves in the vicinity of 
the only joint intersection point that appears to be triv- 
ial: A c = 0. Size scaling of the ratio MPN/N (see Fig. 03) 
reveals no transition too; all MPN/N size-scaling curves 
for non-zero magnitude of disorder decrease with system 
size, as they do for a localized (or critical) state. Thus, 
no signatures of the LDT can be observed in the marginal 
case, indicating that all states are localized. 

In summary, we studied numerically the critical prop- 
erties of the ID two-parameter tight-binding model with 
diagonal disorder and non-random long-range interac- 
tion, J mn = J/\m — n| M , J > and 1 < \i < 3/2. 
The transition point was detected by means of the level 
and wave function statistics. We used the conjecture on 
the scale invariance of the distribution function of the 
nearest-level spacing and the participation number at 
criticality. We find, in particular, that the critical point 
for fi = 4/3 is A c = (11.09±0.21) J. In the marginal case 
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FIG. 5: Size scaling of the ratio MPN/N calculated for fi — 
3/2 and different disorder magnitudes. Thin horizontal dotted 
line is a guide for the eye. 

(/j = 3/2), that is analogous to the standard 2D Ander- 
son modeljSi the only joint intersection point is A c = 0, 
indicating that all states are localized for a finite disor- 
der. 



We demonstrated that finite size effects are very pro- 
nounced within the considered model. Level statistics ap- 
pears to be more affected by these effects as compared to 
the participation number statistics. The dominant con- 
tribution to finite size effects is determined by the size 
dependence of the band width. To obtain the critial dis- 
order, we use a reformulated finite size scaling procedure 
that is corrected for irrelevant size dependencies. 

To conclude, we stress that the scale-invariance of the 
relative fluctuation of the participation number at tran- 
sition is a consequence of critical wave function fluctua- 
tions. We conjecture, therefore, that the analysis of the 
relative fluctuation of the participation number provides 
a general tool to monitor the LDT. The proposed method 
proves to work well for the standard 3D Anderson model 
tooi^i We believe also that this property holds at the 
mobility edge, allowing therefore to monitor the latter. 
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